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The conventional series in powers of the coupling in perturbative QCD have zero radius of conver- 
gence and fail to reproduce the singularity of the QCD correlators like the Adler function at as = 0. 
Using the technique of conformal mapping of the Borel plane, combined with the "softening" of 
the leading singularities, we define a set of new expansion functions that resemble the expanded 
correlator and share the same singularity at zero coupling. Several different conformal mappings 
and different ways of implementing the known nature of the first branch-points of the Adler function 
in the Borel plane are investigated, in both the contour-improved (CI) and fixed-order (FO) versions 
of renormalization-group resummation. We prove the remarkable convergence properties of a set of 
new CI expansions and use them for a determination of the strong coupling from the hadronic r 
decay width. By taking the average upon this set, with a conservative treatment of the errors, we 
obtain as{M^) = 0.31951^1^. 



PACS numbers: 12.38.Bx, 12.38.Cy 

I. INTRODUCTION 

The conventional perturbation expansion of the QCD 
correlators in powers of the strong coupling is prob- 
lematic, because the function that is expanded, like for 
instance the Adler function Z), viewed as a function of 
the coupling, is known to be singular at the expansion 
point, as — [1]. On the other hand, the powers a" 
are holomorphic and, therefore, they can tell us nothing 
about the singularities of the expanded function, includ- 
ing their very existence. As a consequence, no finite-order 
perturbative approximant can share the singularity with 
the expanded function at zero coupling. Singularities can 
emerge only from the infinite series as a whole, which, un- 
fortunately, is not defined, since the perturbation series 
is divergent 0-111 • 

A perturbation series would be more instructive if the 
individual finite-order approximants could retain at least 
some information about the known singularities of the 
QCD correlators. Such approximants would, in every 
order of approximation, bear some information about 
the singularities of the expanded function and, moreover, 
would tell us more about the function from the numerical 
point of view. 

An approach proposed by us some time ago @ con- 
sists in replacing the conventional set of powers of as 
(occurring in the standard perturbation expansion) by a 
new set of functions, in which the available information 
about the singularities of the expanded function is built 
in. In order to define such a new perturbation series, 
two methods can be used: (A) extension of the region of 
convergence by conformally mapping the region of holo- 
morphy of the expanded function onto a disk and 
(B) singularity softening, discussed for instance in 
When combined in a suitable way, they are mutually am- 
plifying their effect. 

The method of conformal mapping was introduced and 
applied in particle physics in [7|, with an intent of ex- 



tending the convergence region beyond the circle of con- 
vergence of an expansion and increasing the convergence 
rate at points lying inside the circle. In the context of 
perturbative QCD, the properties of the new expansions 
based on this technique were investigated in [y, |^, ITol |. 
while in [llj the method was applied for a determination 
of tts from hadronic t decays. 

The conformal mapping method Q is not applicable 
to the (formal) perturbative series of D in powers of as , 
because D is singular at the point of expansion^. On the 
other hand, the method can be applied @, rather than 
to Z?, to its Borel transform B{u) (the precise definition 
of this function will be given in Sec. |lTl. Being holomor- 
phic in a region containing the origin u = of the Borel 
complex plane, B{u) can be expanded in powers of the 
Borel variable as 

n>0 

This series is convergent inside the circle centered at the 
origin u — Q and reaching the nearest singularity of B[u). 
It often happens in practice that the disk of convergence 
of (m is considerably smaller than S, the holoniorphy 
domain of B{u) (we assume throughout the paper that 
B is simply-connected). It is the method of conformal 
mapping (A) that can significantly extend the region of 
convergence of (IT|) , by replacing the series in powers of u 
with a scries in powers of a new variable, w{u)^ 

B(zi) = ^c„(i5(ii))". (2) 

n>0 



^ In the so-called "order-dependent" conformal mappings, which 
were defined also in the coupling plane [12| . , the singularity is 
shifted away from the origin by a certain amount at each finite- 
order, and tends to the origin only when an infinite number of 
terms are considered. 
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Here w{u) denotes the function that conformaUy maps 
the holomorphy region B in the u-plane onto the unit 
disk \w\ < 1 in the w = M;(u)-plane centered at the ori- 
gin (the explicit form of w(u) will be given in Sec. IIIip . 
The expansion coefficients c„ are determined by the 6„ 
coefficients and by the conformal mapping function w(u). 
The expansion ([2]), unlike ((IJ, is convergent in the whole 
region of holomorphy B. Moreover, as proved in Ref. [Tj 
and will be discussed also in Sec. |nT]of the present paper, 
it provides the fastest large-order convergence rate. 

In the present work we focus on the procedure of sin- 
gularity softening, which exploits the known nature of 
the leading singularities of the correlators in the Borel 
plane, by coinpensating them with suitable factors. As 
discussed in 0, [ll|, , this procedure is not unique, in 
contrast with the definition of the optimal mapping jTj , 
which is unambiguous. 

To illustrate the method, we consider the Adler func- 
tion in massless QCD. In the next section we briefly re- 
view the properties of this function in perturbative QCD. 
In Sec. mil we discuss the method of conformal mapping 
for improving the convergence rate of power series 7]. 
Although the method has been adopted and applied by 
many authors, its mathematical foundation is not very 
widely known. We therefore formulate two lemmas which 
allow the definition of the so-called " optimal" conformal 
mapping in this context. The proof of the lemmas is 
presented in Appendix. 

In Sec. llVi we discuss various possibilities of exploit- 
ing the known nature of the leading singularities of the 
Borel transform, and in Sec. |Vl we define a class of new 
expansion functions that implement the two ingredients, 
singularity softening and expansion in a new variable. 
Both the contour-improved and the fixed-order versions 
of the expansions are presented. In Sec. I VII we illus- 
trate the convergence properties of the new expansion 
functions using some mathematical models for the Adler 
function. 

The determination of the strong coupling as is one of 
the most important tests of QCD. As discussed in [l5| . 
the recent determinations at various scales are in an im- 
pressive agreement among each other. The hadronic de- 
cays of the r lepton provide one of the most precise deter- 
mination, which is particularly interesting as it concerns 
a relatively low scale, the mass of the r. The recent cal- 
culation of the Adler function to four loops j l6 [ . the same 
order to which the /3- function is known (iTI . llSll. renewed 
the interest in the determination of aJM3) [l9j|-[28l|. In 



Sec I VIII we present an updated calculation of the strong 
coupling using the new CI expansion functions defined in 
this work. Finally, Sec. I Villi contains a summary of the 
work and our conclusions. 



II. ADLER FUNCTION 

We consider the Adler function [l^, i.e. the logarith- 
mic derivative of the correlation function of two hadronic 



currents, which is expressed in massless perturbative 
QCD as 

D{s) = ^[if„ + k„(-s/m2)] (a,(M2))", (3) 



where s = is the momentum variable and as{fi^) = 
as{p'^)/TT is the strong coupling at the renormalization 
scale /i^. From studies of classes of Feynman diagrams 
it is known that the series in the right hand side 

of ^ is divergent. This series is often assumed [1| to be 
an asymptotic expansion in the limit — >■ 0. Then the 
equality sign in (jS]) is interpreted as the sign used for 
asymptotic expansions. 

The first coefficients Kn calculated in the MS scheme 
are 

Ki = 1, K2 = 1.63982, K3 = 6.37101, = 49.0757. 

(4) 

Several estimates of the next coefficient are avail- 
able: the value Kr, = 378, obtained from the assump- 
tion of a geometrical growth, was adopted in [l^ , while 
the "Fastest Apparent Convergence" (FAC) principle [3^ 
predicts = 275 [l^, |23 [■ A slightly different value, 
— 283, was adopted in 2l|. 
The coefficients Kn{—s/ii^) depend on the 
renormalization-group (RG) /3-function, which is 
calculated at present to four loops [TtI . |Tg. The first 
coefficients /3j in the MS scheme for n/ = 3 are 

;3o = 9/4, /3i= 4, /32 = 10.0599, /^s = 47.228. (5) 

An additional term is sometimes added to the perturba- 
tive expansion of the /3-function, assuming a geometrical 

growth, /34 = ±/3|//32 'mm- 

The Adler function plays a crucial role in the determi- 
nation of as(M^) from hadronic r decays. The method 
is discussed in the seminal paper 31 1 and is reviewed in 
several recent articles [l^ [2l|, [2a, [27 1 . For completeness 
we give below a few details. 

The inclusive character of the total t hadronic width 
makes possible an accurate calculation of the ratio = 
T[t~ — ^ ^'T-hadrons]/r[r^ — > ^-T-e^Fe]. Of interest is 
the Cabbibo allowed component which proceeds either 
through a vector or an axial vector current, since in this 
case the power corrections are especially suppressed. It 
can be expressed in the form 



Rt. 



V+A 



EW 



\Vud\ 



1 + -^f") + <5U + 5- 
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(6) 



where Nc = 3 is the number of colors, S'ew and S'-^^ 
are electroweak corrections, 5pc denotes nonperturba- 
tive power corrections which arise in the framework of 
the operator product expansion (OPE), and ^^"^ is the 
genuine perturbative QCD correction. Unitarity implies 
that this quantity can be written as an integral over the 
spectral function of the polarization function along the 
timelike axis. As shown in 3lJ, the analytic properties of 
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the polarization function and Cauchy theorem ahow one 
to write equivalently (5^"-* as the contour integral 



2m J s \ M2 J 

\s\=M^ 



1 



in) ^<"'' 



(7) 

As discussed in [31|, perturbative QCD is valid outside 
the timelike axis, so the Adler function can be calculated 
along the complex contour using the expansion ([3]). 

More generally, of interest are the moments of the spec- 
tral function, defined for arbitrary sq either as (32j 



or as l2C 



1 



ds 
s 



1 - 



k+l 



D{s), k>0, (8) 



1 

2TTi 



ds 
s 



- 1 



fc > 1. 



(9) 



s|=so 



The main ambiguity in the evaluation of these con- 
tour integrals is related to the renormalization scale. The 
choice ^2 = — s, when ([3]) reads 



5(s) = ^X„(a,(-s))", 



(10) 



n>l 



leads to the so-called " contour- improved" (CI) expan- 
sion [s?! , where the coupling is determined by solving the 
renormalization group equation exactly along the circle. 
The more conventional fixed-order (FO) expansion Q of 
D{s), when /i^ = (or more generally /i^ = sq), is 
obtained formally from (|10p by expanding the running 
coupling as{—s) as 



(11) 



The coefficients depend on the parameters Pk, k < j, 
and the powers of ln(— s/M^), which can acquire large 
imaginary parts for s on the integration circle near the 
timelike axis. As discussed recently 0,111, HI], the dis- 
crepancy between the results given by the CI and FO 
expansions is the main theoretical error in the extraction 
of a,(M2). 

As already mentioned, the renormalized perturbation 
series ^ or ([TU)) are divergent, the coefficients displaying 
at large orders a factorial growth, if„ ~ n\. From inde- 
pendent arguments it is known that correlation functions 
like D, regarded as functions of as, are singular at = 
[l|. For QED, where these facts are well-known [sj], the 
divergence of the series does not affect the phenomeno- 
logical predictions since the coupling is very small. By 
contrast, for a large coupling like q;<,(M^) in QCD the 
consequences are nontrivial. 

The information about the high-order behavior of the 
series ^ is included in the singularities of the Borel 



transform B{u), defined by the series H]), with the co- 
efficients bn related to Kn by 



+1 



n > 0, 



(12) 



where /3o is the first coefficient of the /3-function given 
in ([5]). According to present knowledge Q, B{u) has 
singularities on the real axis for u < — 1 and u > 2, 
known as ultraviolet (UV) and infrared (IR) renornialons, 
respectively. In the present paper we assume there are 
no other singularities in the complex plane, so that the 
holomorphy domain B is the it-plane cut along the real 
axis for u < —1 and u >2. 

The expansions ^ and ([T0| can be formally obtained 
from B{u) by means of an integral of Borel-Laplace type. 
The recovery of the function D(s) is actually ambiguous: 
there are many integral representations admitting ([3]) or 
(|10p as asymptotic expansions (for a recent discussion, 
see [3^). As shown in 3.0], the definition based on the 
principal value prescription. 



D{s) 



1 



■PV 



e /3o»=(-=>) B{u) du, 



(13) 



yields a function D(s) satisfying to a large extent the 
general analyticity requirements in the s-plane, and we 
shall adopt this definition. 



III. ACCELERATING CONVERGENCE BY 
CONFORMAL MAPPINGS 

Because of the first UV singularity at u = — 1, the ex- 
pansion ([1]) converges only on the disk |u| < 1, although 
B{u) is holomorphic in a much larger region. The domain 
of convergence and the convergence rate can be increased 
by expanding the function in powers of a different vari- 
able, defined by the conformal mapping of B (the cut 
M-plane), or a part of it, onto a disk (without loss of gen- 
erality the disk can be taken of radius equal to unity, and 
we shall adopt this convention). It may intuitively seem 
that the larger the domain mapped onto the unit disk, the 
better the convergence properties of the series expansion 
in powers of the new variable. This is indeed true, and 
we shall give this hope a precise mathematical form. The 
result, proved in Ref. [t'I, important and interesting as it 
is for a number of applications, did not raise enough in- 
terest as it deserved, in spite of the many applications of 
the conformal mapping method during the decades. We 
shall therefore state below the main ideas of the proof, 
in order to make the present paper self-contained. The 
following two lemmas show which is the variable that 
provides the best asymptotic rate of convergence. 

Lemma 1: Let 2?i and I?2 be two domains in the complex 
M-plane, with I?2 C 7^ 2^1, such that the confor- 

mal mappings 2?i — > ICi and 2?2 ~^ ^2 exist, where ICi 
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and IC2 are unit disks. Consider the two conformal map- 
pings 

zi = zi{u) : 2?i ^ /Ci = {zi : Izil < 1}, 

Z2 = hiu) : 2?2 ^ /C2 = {Z2 : \z2\ < 1}. (14) 

Let Q be a point of I?2, Q G I?2, such that zi{Q) = 
and hiQ) = 0. Then 

|zi(u)| < |z2(m)|, forah uel?2, u^^Q. (15) 

Lemma 2: Let I?i and I?2 be the domains defined in 
Lemma 1, zi(u) and 22 (w) the mappings (|14p and B{u) 
a function holomorphic in 2?i. Define the expansions 



= X]c„,i(zi(^i))", 



n=0 



B(7X) =^C„,2(52(w))", 



(16) 



(17) 



which are convergent for zi = zi{u) G /Ci and Z2 = 
Z2{u) € JC2, respectively. Assume in addition that the 
hmits hm„^oo \/|c„,i| and hm„^oc. \/\cn,2\ exist^ and 
are equal to one: 



lim \ \cn.\\ = lim \/|c„,2| 

n^oo V n— )-oo V 



1. 



(18) 



Then a positive integer N = N(u) exists such that the 
following inequality holds: 



c„4(zi(u))" 



c„,2(22(u))" 



< 1, 



(19) 



for any n integer, n > N, and m € I?2, u ^ Q. 

Proofs of Lemma 1 and Lemma 2 are presented in Ap- 
pendix. 



boundary of B is mapped on the boundary circle of the 
unit disk and the series ^ in powers of w{u) is conver- 
gent everywhere in B. Moreover, the asymptotic con- 
vergence rate of this series is, at any point u £ B, the 
fastest among all conformal mappings. This mapping is 
known as "optimal" conformal mapping for convergence 
acceleration 0, [HI . 

Let us discuss an example to illustrate this result. As- 
sume we decide to modify this mapping by adding a re- 
gion lying outside the holomorphy domain. In doing so, 
we wedge a region containing singularities inside the cir- 
cle and unavoidably make the convergence radius smaller. 
As a consequence, the large-order convergence rate is 
worse. 

If, on the other hand, we omit to map a part of the 
holomorphy region inside the unit circle, the convergence 
rate becomes worse, as follows by a direct application of 
(fT91) . The reason is that we do not make full use of analyt- 
icity in this case, leaving aside a part of the holomorphy 
region. We conclude that by w(u), the optimal confor- 
mal mapping function, neither any singularity is mapped 
inside the circle, nor any part of the holomorphy region 
is left out of the circle. 

As can be seen from the proof of Lemma 2, when n is 
large enough, the inequality (|A8() reduces to \np{u) < 0, 
i.e. the coefficients c„j play no role in the ratio TZn{u) 
of the convergence rates. On the other hand, when n is 
finite (and small), the term g{n)/n in (|A8[) may be pos- 
itive and greater than In p{u) , making the ratio TZn (u) 
greater than one. This may happen, in particular, when 
both |zi(u)| and |22('«)|, as well as their ratio p, are close 
to 1. Therefore, the expansions in powers of other con- 
formal mappings can provide at finite orders a better ap- 
proximation compared to the optimal expansion, a fact 
observed numerically in some cases, especially for points 
near the boundary of the analyticity domain. 

For the Adler function, assuming that there are no 
other singularities except for the cuts along the real axis 
for M < — 1 and u > 2, the optimal conformal mapping is 

i 



A. Optimal conformal mapping 

It is now easy to understand from Lemma 1 and 
Lemma 2 that the larger the domain mapped onto the 
disk, the larger the domain where the expanded function 
is represented by a convergent power series and, also, 
the faster the convergence rate at a given point. From 
the inequality p^ . it follows that the best asymptotic 
convergence rate is obtained with the variable {^{u) that 
maps B, the holomorphy domain of B(u), onto the unit 
disk |ui| < 1 in the plane w = {^{u). In this case, the 



w{u) 



with the inverse 



u{w) 



8w 



3~2w + 3w^' 



(20) 



(21) 



Using the optimal expansion ^ and the definition 
(jl3L we were led in a natural way to the new pertur- 
bative expansion Q 



n>0 



(22) 



^ The essence of this is that the expansions 1II6II and II17II have 
equal radii of convergence. This assumption is nontrivial, be- 
cause the expanded (Adler) function might be of such a form 
that certain singularities of B{u) in 21 or 22 might disappear. 



W„(s) = -PV 
Po 



>o»=(-=)(w(u))"du. (23) 
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By construction, the series when reexpanded in 

powers of as, reproduces the expansion (jlOp with the 
coefficients Kn known from Feynman diagrams. On the 
other hand, the expansion functions W„ are singular at 
tts = 0, resembhng the expanded function D itself [l0| . 
Moreover, as discussed in @, under certain conditions, 
the expansion (j22p converges in a domain of the s-plane. 



IV. SINGULARITY SOFTENING 

In the particular case of the Adler function in massless 
QCD, the nature of the leading singularities in the Borel 
plane is known: near the first branch points, u — —1 and 
u = 2, B{u) behaves like 



Biu) 



(1 + up^ 



and B{u) 



r2 



(l-u/2)'r2' 



(24) 



respectively. The residues ri and r2 are not 
known, but the exponents 71 and 72, calculated us- 
ing renormalization-group invariance, have known values 



71 



1.21, 



72 = 2.58 . 



(25) 



The expansion ^ takes into account only the position of 
the renormalons in the Borel plane. If a sufficient num- 
ber of expansion coefficients were known, ([2]) would be 
expected to describe also the character, strength, etc., 
of the singularities as well. Since, however, only a few 
perturbative coefficients are at present explicitly avail- 
able, one cannot expect that the expansion of the type 
might be able to give a satisfactory approximation of 
B{u). It is better than (jlj, which has no singularities in 
any finite-order approximation, while every finite-order 
approximant to ([1]) has the same location of cuts as the 
expanded function. But it can hardly be expected that 
the first four or five perturbative coefficients would be 
able to represent B{u) with a satisfactory accuracy. 

An explicit account for the leading singularities ((24)) 
would therefore be helpful to further improve the con- 
vergence. This can be done by multiplying B{u) with 
suitable factors that vanish at u = — 1 and u — 2 and 
compensate the dominant singularities. The subsequent 
expansion of the product in powers of a conformal map- 
ping variable is expected to converge better. This proce- 
dure is known as " singularity softening" 0, Q . 

In contrast with the optimal conformal mapping, sin- 
gularity softening is not unique [ll|, [U^ HI] • The singu- 
larities are present in B{u), but we do not know their 
actual form, except for the behavior (1^^ near the cor- 
responding branch-points. A possibility is to multiply 
Bhi) by simple factors like (1 + u)i^ (1 - u/2)'<^ In 
[Tl| the alternative softening factors (H- u;)^'''^ (1 — w)^'^^ 
were adopted, where w = w{u) is the optimal mapping 
(pn|. The product of B{u) with these factors was after- 
wards expanded in powers of the same variable w. Other 
possibilities will be investigated in the next section. 



V. NEW EXPANSION FUNCTIONS 

The product of B{u) with softening factors is expected 
to contain milder singularities, which vanish instead of 
exploding at u = — 1 and u = 2 (in very peculiar cases 
the singularities may disappear altogether, but this sit- 
uation is very unlikely). The effect of a mild singular- 
ity in a function is not visible at low orders in its series 
expansions, and is expected to appear only at large or- 
ders. Therefore, we can ignore their effects, expanding 
the product in powers of variables that account only for 
the next branch-points of B{u). In the case of the Adler 
function, these singularities are placed at it = 3, 4, etc., 
on the positive axis, and at m = —2, —3, etc., on the 
negative axis. 

In general, we consider the functions 



which map the w-plane cut along u < —j and m > fc to 
the disk \wjk\ < 1 in the plane wjk = Wjk{u). For j = 1, 
fc = 2, we recover the optimal mapping (j20p . In the fol- 
lowing, we shall consider also the variables W12, W13, wioo 
and 1023, for which the corresponding unit disks \wjk\ < 1 
are shown in Fig. [H We mention that the mapping wioa, 
suggested in [2|, was discussed in a similar context in 
[39| . and W13 was investigated in According to the 
discussion in the previous section, the last three map- 
pings "push" inside the unit circle a part of the u-plane 
containing some singularities (indicated in Fig. [T]). As 
a consequence, the expansions based on these variables 
will converge in a smaller domain and their convergence 
rates will be, in principle, worse than that of the optimal 
mapping ■Wi2. 

According to the above discussion, we shall expand in 
powers of Wjk the product of B{u) with suitable softening 
factors. Specifically, we consider the expansions 



(27) 



where Sjk{u) must "soften" in principle all the singular- 
ities of B{u) at — j < u < and < u < k. Numerically, 
it appears to be convenient to choose the factor Sjk as a 
simple expression with a rapidly converging expansion in 
powers of Wjk , thus ensuring a good convergence of the 
product ([77)1 . 

A systematic application of this idea to the singulari- 
ties of B{u) requires the knowledge of the nature of the 
branch-points, which at present is limited to the lead- 
ing singularities at u = — 1 and u — 2. Therefore, we 
shall limit ourselves to compensating factors that vanish 
at these points, and take Sjk of the form: 



Sjkiu) = 1 - 



Wjfe(-l) 



1 - 



Wjkju) 

Wjk{2) 



(28) 
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Re w,. 




FIG. 1: The unit disks \wjk\ < 1 on which the conformal mapping defined in (|26|l maps the cut it-plane, for several values 
of j and k. In the last three figures, the thick lines indicate the residual cuts along the segments (wi3(2), 1), (wicx)(2), 1) and 
(— 1, W23(— 1)), (u)23(2), 1), which correspond, in the xt-plane, to the segments (2,3), (2,cx3), and (—2, —1), (2,3), respectively. 



The exponents 7^ = 71(1 + and ^2'^' =72(1 + 42), 
where 6ij is Kronecker's function, are taken such as to 
reproduce the nature of the first branch-points of B{u), 
given in ((M)) . In particular, for the optimal case j = 1, 
k = 2we recover from ^ the factor (1 + w)2ti(1-w)272 
used in [iJl, with w — w{u) defined in (pn| . 

Strictly speaking, for a fixed pair (j, k) the ex- 
pansion (j27p converges only on the disk \wjk\ < 
min[|-u;jfc(— 1)1, |wjfc(2)|]. For the optimal choice j — 
1, fc = 2, the expansion converges in the whole unit disk 
1^12! < 1, i.e. in the whole w-plane except for the cuts 
along the real axes for u > 2 and u < — 1 Q. For other 
mappings, the convergence disk is limited by the begin- 
ning of the cuts shown in Fig. [1] In particular, if j = 1 
and fc > 2 the expansions (j27p diverge for real u greater 
than 2, while for the conformal mappings with j > 1, the 
expansions start to diverge for u greater than one, due 
to the singularity at m = — 1 pushed inside the circle (as 
in the last case shown in Fig. [T]). However, the prod- 
uct Sjk{u)B{u) has only mild singularities. Moreover, 
the expansion (P7)) enters the Laplace-Borel integral 
where, for values of in the domain of interest, the con- 
tribution of high values of u is suppressed. In particular, 
if as is not very large, the region u > 2 brings a small 
contribution to the integral, so signs of divergence in the 
case of the variables W13 and wioo are expected to occur 
only at very large orders N . On the other hand, for the 
variable W23, it is natural to expect signs of divergence at 
lower values of N , since the series (P7)) does not converge 
for u > 1 . 

By combining the expansion (P7)) with the definition 
(ITSl) . we are led to the class of expansions 



(k) 



The new fixed-order (FO) expansions can be obtained 
in a straightforward way [l Ij , using as starting point ([3|) . 
They have the generic form 



ra>0 



Sjk{u) 



du. 



(31) 



(32) 



The expansion functions W}^ coincide with the opti- 
mal functions investigated in detail in [ll|. In the fol- 
lowing, we shall also consider the expansion in terms of 
the functions W}^ , and (and their correspond- 

ing FO versions). We emphasize that these expansions 
contain different softening factors, which coincide only 
for u near the corresponding singularities, when they re- 
produce the known behavior ((24| . The treatment of the 
residual singularities after softening is also different: the 
expansion in powers of the optimal mapping uses the po- 
sition of the first singularities, which in general do not 
disappear completely after the multiplication with the 
compensating factors. The other expansions exploit the 
fact that a mild singularity can be neglected at low per- 
turbative orders, and use also some information about 
the position of higher renormalons. So, the representa- 
tions for different j and fc can be considered 
independent perturbative expansions of the Adler func- 
tion. 



D{s) 



n>0 



Wr(s) = -PV /e-^w-^ 



S]k{u) 



(29) 



(30) 



By inserting into ([50]) the coupling as{—s) calculated by 
solving the renormalization-group equation for s along 
the circle defined in the integral ([T]), we obtain the 
countour-improved (CI) version of the new expansions. 



VI. MODELS 

For testing the convergence of the various expansions, 
we consider a class of models of the type proposed in 
[21I I ((but analyzed without using conformal mappings), 
which parametrize the Borel transform B{u) and then re- 
cover the Adler function by means of the PV prescription 

In the model proposed in [21j, the function B{u) is 
expressed in terms of a few UV and IR renormalons 

Bbj(u) = BUV(^)+BiR(^)+BiR(„)+dPO+rfPO^^ (33) 
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FIG. 2: Real part of the Adler function of the model [SJ defined in (|33p -(l36 ( l. calculated along the circle s = exp{i(l>) 
for Qs(M^) = 0.3156, using the perturbative expansions with A*" = 5 terms. Left panel: CI expansions. Right panel: FO 
expansions. The exact function is represented by the solid line. 
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FIG. 4: Real part of the Adler function of the alternative model defined in (|37p -(l39 p . calculated along the circle s — exp(i(j!)) 
for asiJsA^) = 0.3156, using the perturbative expansions with A = 5 terms. Left panel: CI expansions. Right panel: FO 
expansions. The exact function is represented by the solid line. 



parametrized as 



7p 



l + fel(p-M)+fe2(p-w)^ + 



(34) 



l + hi{p + u) + b2{p + uf + .. 



(35) 

Most of the parameters are fixed using a renormalization- 
group analysis at four loops, the free parameters of the 
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FIG. 5: As in Fig. |4]for Af = 18. The standard CI and FO expansions exhibit big oscillations and are not shown. 



TABLE I: The quantity 5'°' for the model Bbj proposed in [Ul and specified in (ES-lIMl), calculated for as(M^) = 0.34 with 
the standard and modified CI and FO expansions truncated at the order A'^. Exact value 5'"' = 0.2371. 

N CI St. FO St. CI wi2 FO wi2 CI wxz FO luia CI wi^o FO luioc CI w^z FO wn 



2 


0.1776 


0.1692 


0.1977 





2228 





2070 





2203 





1883 





2524 





2123 





2099 


3 


0.1898 


0.2026 


0.2009 





2460 





2030 





2440 





1975 





2530 





2028 





2437 


4 


0.1983 


0.2200 


0.2263 





2463 





2194 





2460 





2288 





2465 





2206 





2463 


5 


0.2022 


0.2288 


0.2290 





2440 





2268 





2423 





2310 





2427 





2292 





2423 


6 


0.2046 


0.2328 


0.2324 





2484 





2306 





2421 





2321 





2431 





2319 





2449 


7 


0.2046 


0.2342 


0.2339 





2536 





2331 





2457 





2333 





2454 





2345 





2502 


8 


0.2017 


0.2353 


0.2339 





2505 





2343 





2484 





2341 





2471 





2347 





2476 


9 


0.2004 


0.2367 


0.2341 





2431 





2348 





2457 





2346 





2465 





2347 





2377 


10 


0.1842 


0.2390 


0.2351 





2420 





2348 





2394 





2348 





2436 





2353 





2337 


11 


0.1962 


0.2402 


0.2359 





2406 





2348 





2352 





2349 





2399 





2348 





2335 


12 


0.1123 


0.2436 


0.2362 





2298 





2351 





2349 





2349 





2370 





2374 





2262 


13 


0.2629 


0.2408 


0.2362 





2229 





2355 





2341 





2349 





2356 





2348 





2226 


14 


-0.2915 


0.2575 


0.2364 





2242 





2361 





2303 





2349 





2354 





2395 





2314 


15 


1.1011 


0.2170 


0.2367 





2173 





2366 





2277 





2350 





2357 





2356 





2365 


16 


-3.362 


0.3818 


0.2368 





2102 





2369 





2305 





2351 





2360 





2343 





2374 


17 


9.5931 


-0.1881 


0.2368 





2176 





2372 





2356 





2352 





2360 





2533 





2512 


18 


-31.52 


2.144 


0.2368 





2201 





2373 





2371 





2354 





2359 





1926 





2665 



models being the residues d}^^ and of the first 
renormalons and the coeficients <SP , d^'~^ of the polyno- 
mial in (j33p . They were fixed in f21'| by the requirement 
to reproduce the perturbative coefficients A'„ for n < 4 
from (HI) and the estimate — 283, and are: 

d^^ = - 1.56 X 10"^ 4^ = 3.16, 4" = -13.5, 

4° = 0.781, 4° = 7.66 X 10"^ (36) 

We also investigated alternative models, where we im- 
posed a specific residue at ti = 2. In one such exam- 
ple, we kept the same expressions as in [2l| for the first 
three singularities and the same values of the residues 
at u = — 1 and it = 3, while choosing a smaller residue 
at M = 2, d^ — 1. The model must contain then three 
additional free parameters in order to reproduce the first 



five Kn- Specifically, we introduced a quadratic term in 
the polynomial and two additional IR singularities, at 
u = A and u — 5. For convenience, the nature of these 
additional singularities, which is not known, was taken 
to be same as that of the u — 3 singularity. Thus, we 
considered the alternative model: 

+ ^^^^+4°+rf^"+4°"^ (37) 

where, as discussed above, we took as input 

^ _ 1.56 X 10-^ 4^ = 1, 4^ = -13.5, (38) 
and determined the remaining five parameters by match- 
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TABLE II: The quantity for the modified model Bait specified in calculated for a^M^) = 0.34 with the standard 

and modified CI and FO expansions truncated at the order A*'. The rows for N < 5 are identical to those in Table [T] Exact 
value S^°^ = 0.2102. 



N 


CI St. 


FO St. 


CI ^12 


FO wi2 


CI Wl3 


FO wi3 


CI Wlao 


FO 


CI li'23 


FO V02d. 


6 


0.2041 


0.2318 


0.2263 


n 




0.2271 


0.2420 


n 


2284 


0.2431 


0.2260 


0.2454 


7 


0.2041 


2290 


0.2201 


n 


2628 


2220 


0.2481 


n 


2230 


0.2472 


0.2174 


2580 


8 


202"^ 


0.2213 


n 2202 


n 


97^R 
z t uu 


0.2164 


2595 


n 

u 


91 89 


0.2541 


0.2136 


0.2734 


9 


0.2037 


0.2110 


0.2175 


n 


2742 


0.2143 


2686 


n 


2154 


2608 


0.2138 


0.2706 


10 


0.1924 


n 2032 


2055 





2709 


0.2144 


0.2651 





2146 


2629 


0.2115 


0.2517 


11 


0.2124 


0.2004 


0.1982 





2905 


0.2136 


0.2504 





2146 


0.2578 


0.2068 


0.2531 


12 


0.1412 


0.2071 


0.2007 





3063 


0.2111 


0.2406 





2148 


0.2468 


0.2081 


0.2627 


13 


0.3121 


0.2117 


0.2022 





2820 


0.2086 


0.2449 





2149 


0.2340 


0.2060 


0.2133 


14 


-0.2105 


0.2344 


0.2001 





2666 


0.2074 


0.2459 





2146 


0.2239 


0.2124 


0.1338 


15 


1.2336 


0.1934 


0.2009 





2865 


0.2079 


0.2176 





2142 


0.2187 


0.2087 


0.1192 


16 


-3.147 


0.3500 


0.2044 





2562 


0.2091 


0.1676 





2136 


0.2175 


0.2073 


0.0930 


17 


9.948 


-0.2333 


0.2059 





1822 


0.2102 


0.1355 





2130 


0.2175 


0.2275 


-0.0415 


18 


-30.94 


2.084 


0.2058 





1722 


0.2107 


0.1345 





2124 


0.2159 


0.1617 


-0.1019 
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FIG. 6: Moment M5 defined in ((8| for the model [IJ], calculated for |so| = and as(M^) = 0.34, as a function of the 
perturbative order N , for the standard and the new expansions based on several conformal mappings Wjk defined in (|26l) . The 
grey horizontal line is the exact value. Left panel: CI expansions. Right panel: FO expansions. 



ing the same coefRcients Kn for n < 5: 

<i^° = 3.2461, = 1.3680, 4° = 0.2785, 

4^ = 1560.614, = -1985.73. (39) 

We emphasize that we consider these models only as a 
mathematical frame to test the convergence properties of 
the various expansions. The physical plausibility of one 
model or another [21, 24] will not be discussed here. 

In Figs. [Hand [21 we show the real part of the Adler 
function for the model 21], calculated along the circle 
s — AI^ exp(i(/)) with the standard and the new CI and 
FO expansions defined in (|29l) - (|32l) . where the perturba- 
tive expansions were truncated at iV = 5 and N = 18, 



respectivel y. T o facilitate the comparison with previous 
works we took as{M^) = 0.3156 in this calcu- 

lation. For TV = 18, the standard expansions exhibit big 
oscillations and are not shown. 

The curves show that the new CI expansions based on 
conformal mappings give a good approximation, which 
improves with increasing TV, of the real part of D{s) along 
the whole circle (only the mapping W23 shows signs of 
divergence for large TV, as expected). This behavior is 
valid also at higher TV (we explored values up to TV = 
25), for both the real and imaginary parts of the Adler 
function. 

As concerns the FO expansions, the description they 
provide is quite good for points close to the spacelike axis. 



10 



0.8 



0.6 



0.2'^^ 







CI standard 








♦ - 


--♦ 


CIw^3 


•- 


-• 


CIWj^ 


T- 


T 


CIW^3 



6 8 10 12 
Perturbative order 



14 



16 



18 



0.8 



0.6 







FO standard 






FOw,, 


o- 


--0 


FOw„ 


o 


-0 


FO W|_ 




- V 


FOw,3 




6 8 10 12 
Perturbative order 



14 



16 



18 



FIG. 7: Moment Ms defined in © for the model [U], calculated for |so| = and OsiM^) = 0.34, as a function of the 
perturbative order A'^, for the standard and the new expansions based on several conformal mappings. The grey horizontal line 
is the exact value. Left panel: CI expansions. Right panel: FO expansions. 




FIG. 8: Moment M5 defined in ([SJ for the alternative model calculated for |so| = and Qs(M^) = 0.34, as a 

function of the perturbative order for the standard and the new expansions. The grey horizontal line is the exact value. Left 
panel: CI expansions. Right panel: FO expansions. 
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FIG. 9: Moment M5 defined in (O for the alternative model (ETJ-dMl), calculated for jsol = and as{M^) = 0.34, as a 
function of the perturbative order for the standard and the new expansions. The grey horizontal line is the exact value. Left 
panel: CI expansions. Right panel: FO expansions. 



(j) — TT, but gradually deteriorates near the timelike axis, ^ = 0. To understand this behavior, we remark that the 
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coupling as{—s) is calculated along the circle as the exact 
solution of the RG equation in terms of as{M^), both in 
the "true" function (IT3t and the CI expansion functions 
([50]) (recall that Ug = ag/Tr). Therefore, the improve- 
ment of the series achieved by the conformal mappings 
is clearly seen along the whole circle in the case of the 
CI expansions. On the other hand, the FO expansions 
are obtained by expanding as{—s) in powers of as(Af^), 
according to ((TT|) . As remarked in [3j], this expansion 
has a poor convergence near the timelike axis, due to the 
appearance of large imaginary logarithms in the coeffi- 
cients. The curves in the right panels of Figs. and 
[3] show clearly the effect of the weak convergence of the 
additional series pT|) involved in the definition of the FO 
expansions. The detailed behavior depends on the con- 
formal mapping: for instance, the expansions based on 
the variables wi^ and wioo provide, for increasing N, a 
good approximation up to points rather close to = 0, 
as shown in the right panel of Fig. [3] 

In Figs. 13] and [5] we repeat the analysis for the alterna- 
tive model (I571) - ([M)) . The perturbative curves in Figs. 0] 
coincide with those in Fig. [2j since the first five perturba- 
tive coefficients of the two models coincide. On the other 
hand, the "true" function is slightly fiatter in the sec- 
ond model, and is better approximated by the standard 
CI expansion with N — 5 terms than it was the model 
shown in Fig. [51 By increasing TV, the new CI expan- 
sions based on conformal mappings again converge nicely 
towards the true function, as is seen in the left panel of 
Fig. [5] The new FO expansions give a good approxima- 
tion near the spacelike axis, but a poor description, even 
worse than for the previous model, near the unitarity cut. 
The conclusion is that, for both models, the new CI ex- 
pansions give a precise approximation of D{s) along the 
whole circle, while the FO expansions give a description 
that deteriorates near the timelike axis. 

For the determination of as{M^), the quantity of in- 
terest is the integral defined in 0. In Table HI we 
give the values of 5'^"^ for the model [21] j calculated with 
the standard and the modified CI and FO expansions, as 
a function of the perturbative order N. To facilitate the 
comparison with similar results reported in [ill [2l| we 
took as{M^) = 0.34. As discussed in [U, at low iV, the 
standard CI expansion gives values systematically lower 
than the true result, while the standard FO expansion 
give s a better approximation. As remarked already in 
[11|, for the same low values of N, the difference be- 
tween the results of the new CI and FO expansions is 
smaller than that of the standard ones. Moreover, as 
seen from Table HI the difference decreases when pass- 
ing from A'^ = 4 to = 5, contrary to what happens 
with the standard expansions. At larger A^, the new CI 
expansions approach the exact value (deviations appear 
only for the expansion based on the conformal mapping 
W23, for the reasons discussed in Sec. IV]). The new FO 
expansions give slightly worse values, however, the map- 
pings Wis and wioo lead to good approximations at large 
N also in the FO case. 



In Table [TTl we present similar results for the alterna- 
tive model (|57| - ([55)) . By construction, the first five rows 
in Tables U and |lT] are the same (but the " true" value is 
now different, (5^°) = 0.2102 instead of 6'-°^ = 0.2371 in 
Table m. The CI expansions based on the mappings Wi2, 
Wi3 and wioo approach at large A^ the exact value also in 
this case. In the FO case, the description is less precise 
and, for the values of A^ considered, only the expansion 
based on the mapping wioo exhibits a good numerical 
convergence. 

Finally, we illustrate the properties of the expansions 
by calculating the perturbative part of the moments of 
the spectral function, defined in ([5]) and Detailed 
studies of the moments have been performed in Refs. 
[20I [26I . I32I . [33I I , especially in connection with the power 
corrections. 

The approximation provided by various expansions re- 
sults from the interplay between the behavior of the se- 
ries and that of the integrand along the circle. Both in- 
tegrands in (|S]) and © vanish on the timelike axis, but, 
while the second integrand suppresses the contribution of 
a region, which increases with k, near the timelike axis, 
the first exhibits oscillations increasing with k along the 
circle, and vanishes also on the spacelike axis for odd k. 
From the behavior shown in Figs. [2][5j we expect there- 
fore a better approximation by the new FO expansions 
of the moments ^ compared to ^ . We recall that the 
new FO expansions give a good description of the Adler 
function near the spacelike axis, but the accuracy deteri- 
orates near the timelike axis, due to the poor convergence 
of the expansion (jlip . For the new CI expansions, a more 
or less comparable description at low orders, depending 
on the specific integrand, and a very good convergence 
at high A^, are foreseen. For the standard expansions, 
the results depend on the fortuitous cancellations of the 
contributions along the circle, as in the case of S^^\ 

This expectation is confirmed in Figs. EEl where we 
show the moments and M5 calculated with the stan- 
dard and the new expansions for the model (|33p - (P^ de- 
fined in [2l[, and for the alternative model defined in 
((57)) - (|55)) . We chose a rather high moment to see clearly 
the difference between the behavior of ([8]) and ([9]) . In all 
cases, we took sq = M^. 

For the first model, the new CI expansions give a very 
good description of both moments, as shown in the left 
panels of Figs. [S][71 The new FO expansions give a rather 
poor description of the moment M5 , but a very good ap- 
proximation of the moment M5 . The right panel of Fig. 
[7] shows that, for an integrand that strongly suppresses 
the region near the timelike axis, the new FO expansions 
provide a very good description. As for the standard ex- 
pansions, at low orders they give a better approximation 
of the moment M5 , for which suitable cancellations of the 
terms along the circle occur. At larger N, both standard 
expansions show large deviations from the true result. 

For the second model. Figs. [SHn] indicate a similar pat- 
tern, with a slightly worse approximation given by the 
new expansions at low orders. For both models, the most 
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impressive feature is the good description of the moments 
by the new CI expansions at large orders (only the expan- 
sion based on the variable W23 exhibits small deviations 
at large N, as expected from the discussion in Sec. N} . 



VII. DETERMINATION OF a,{M?) 

The above analysis demonstrated the good conver- 
gence properties of the new contour-improved (CI) per- 
turbative series based on singularity softening and expan- 
sions of the Borel transform in powers of suitable confor- 
mal mappings. We apply now these expansions for a 
determination of a^iM'^) from the experimental rate of 
hadronic r decays. We emphasize that our calculation is 
not based on the models discussed in the previous sec- 
tion, but relies only on the known coefficients Kn given 
in (HI) and a very conservative choice [2l|, [2^ for the next 
coefficient, = 283 ± 283. For the running of the cou- 
pling, we use the calculated perturbative coefficients /3j 
from ([5]), an assumption about the next coefficient /34 
being considered only for the assessment of the errors. 

The standard determination of as{M^) from hadronic 
r-decays requires the theoretical calculation of the inte- 
gral defined in ([T]), using the perturbative expansion of 
the Adler function. On the other hand, the quantity (5'^*'^ 
can be determined with great precision from ^ . The re- 
cent determination Rt.v+a = 3.4771 ± 0.0084 El leads 



to the updated phenomenological value [28 1 



phen 



0.2037 ± 0.0040cxp ± 0.0037pc, (40) 



where the first error is experimental and the second ac- 
counts for the power corrections. 

Using this input, the values of as(Af^) obtained with 
the new CI expansions defined in (P^ -([5(7 | . with the 
expansion functions W^^, W^^, W^°° and W^^, respec- 
tively, are: 



0.3195 ± 0.0034exp ± 0.0031pc to;o?37(if5) 
0.3208 ± 0.0035oxp ± 0.0032pc 



+0.0131/7^ N 



+0.0018^ ,N 
-0.0019 v^'-"'^^/' J 

+0.0024/ J N 
-0.00881'''''^^^/' 



0.3182 ±0.0033exp± 0.0031pc to:om(^5) ^^^^(scale), 
0.3193 ±0.0034exp±0.0031pc to onsl-^s) tSilJEJi (scale). 

(41) 

The first two errors are produced by the uncertainties of 



"phen given in (j40l) . the third one is obtained by varying 
the coefficient in the conservative range mentioned 
above, and the last error accounts for the variation of 
the scale as £,M^, with ^ in the range 0.5 — 1.5 [26l |. 

The largest errors in (j4ip are produced by the uncer- 
tainty in the coefficient K^. To understand this result, 
we remark that the series ([^ - ([501) . when reexpanded in 
powers of Os, generate an infinite number of higher-order 
terms 0, [ll|. In particular, the representations based 
on the expansion functions W^"^, W^^, W^°° and W^^, 
truncated after four terms (i.e. neglecting the fifth term 



proportional to K^), lead to coefficients K5 equal to 256, 
161, 256, and 179, respectively. If a very different value, 
like K5 = or K5 — 566, is imposed, the expansions 
can match the same S'^^^^ only with the price of a much 
larger/smaller coupling, respectively. In fact, if K^, is as- 
sumed to be negative and large, the coupling should be 
so large that the calculation becomes unreliable, and no 
solution as{M^) exists at all. 

By taking the average of the values in (|4T|) . we obtain 

a,(M2) = 0.3195 ± 0.0034exp ± 0.0031pc laon'C^O 
l°:°°i^(scale)± 0.0005^,, (42) 

where we added an uncertainty to account for the trun- 
cation of the /3-function (obtained by including a fur- 
ther term based on a geometrical growth, /34 = ±/3|//32 



[19|, |26|). We emphasize that the errors quoted in (|42)) 
were obtained as simple averages of the individual errors 
given in (|4ip. Much lower uncertainties would be ob- 
tained if standard statistical procedures for independent 
determinations (for instance, Eqs. (14) and (15) of |1,5;]) 
were applied. In practice, although the values given in 
(|4ip may be considered independent theoretical determi- 
nations, we prefer the conservative errors given in (j42L 
which avoid any bias. The remarkable consistency of the 
theoretical determinations (HIT) is nevertheless a strong 
argument in favor of our predictions. 

By combining in quadrature the errors given in (j42p . 
we finally obtain 



n ^i^Q'^ +0.0189 

U.OiJd _o.0138- 



(43) 



The central value in (1431) co incides practically with our 
previous determinatioii[ll| , as{M^) = 0.320 ± 0.011, 
obtained with the optimal mapping W12 and the slightly 
different value S'^^^ = 0.2052 ± 0.0050 from The 
smaller error quoted in [11] is due mainly to a smaller 
range, = 283 ± 142, adopted there for the coefficient 

K5. 

We note that for the same input, the standard CI ex- 
pansion to 5-loops leads to as{M^) = 0.3419 ± 0.012, 
while the standard FO expansion gives as{M^) = 
0.3199lg;!5o7^ [H. The smaller errors are mainly due 
to the fact that the standard expansions are less sensi- 
tive to the variation of K^. However, these expansions 
have the behavior expected for an asymptotic series, ap- 
proaching the expanded function up to a certain order 
N, and starting to oscillate violently afterwards. For 
some expanded functions, the minimal error reached be- 
fore the onset of oscillations may be rather small, but in 
other cases, the standard expansion never describes the 
function with sufficient accuracy. Therefore, the uncer- 
tainty of can generate only a part of the truncation 
error, an additional term being necessary in order to ac- 
count for the divergent pattern (this term may be taken, 
for instance, as the difference of about 0.022 between the 
predictions of the standard CI and FO expansions to 5- 
loops). 
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VIII. SUMMARY AND CONCLUSIONS 

In this paper, we investigated a new class of expansions 
of the Adler function in perturbative QCD and apphed 
them to a determination of as from hadronic r decays. 
Our work extends previous studies reported in [ol-flll 

As remarked in (13| . if a series is divergent and the 
expansion parameter is not very smah, a summation of 
the perturbative expansion is indispensable. The defini- 
tion of the expansion functions investigated in this work 
exploits the information available on the large-order be- 
havior of the perturbative series, together with mathe- 
matical results on accelerating the series convergence by 
conformal mappings Q- These techniques are suitable 
for the Borel plane, where an analyticity domain around 
the origin exists. An important feature [9] is that the ex- 
pansion functions share the singularity of the expanded 
function at the origin ctg = of the coupling plane. 

In the present study, we focused on the procedure of 
"singularity softening," by which the strong leading sin- 
gularities in the Borel plane are turned into milder singu- 
larities, where the function vanishes instead of exploding. 
In practice, this is achieved by expanding the product 
of B{u) with suitable factors that vanish at the points 
M = — 1 and u ~ 2. The procedure is possible due to 
the exact results available on the nature of the leading 
singularities ^21,, .38] . Since the effect of a mild sin- 
gularity is expected to occur only at higher orders in a 
power series, one can choose as expansion variable a con- 
formal mapping that accounts only for the further sin- 
gularities of the Borel transform. Extensive numerical 
studies showed that it is convenient to take the factors 
that multiply the function B{u) (which in principle are 
arbitrary) as simple expressions of the same variable that 
is used for expanding the product. 

In Sec. I VII we investigated in detail the properties of 
the expansions defined in Sec. |V] by using two specific 
models for the Adler function. As already mentioned, 
we consider these models only as a mathematical frame 
for testing the convergence properties of the various ex- 
pansions and make no assumption about their physical 
plausibility (actually, there is no consensus on this sub- 
ject in the literature). In all cases, we obtain a good con- 
vergence of the new contour-improved (CI) expansions 
defined in (pn)) - (|30p . for the choices oi j and k adopted in 
this work. The alternative fixed-order (FO) expansions 
(l5T]) - ([5^ converge near the spacelike axis, but provide a 
worse approximation near the timelike axis, since they 
also involve the expansion (jlip . which converges slowly 
in this region. 

Suitable cancellations between the two terms in the 
expansion (|3]) make the standard FO expansions more 
suitable for calculating integrals like ([7]) for some mod- 
els. Such fortuitous cancellations are not expected to 
occur in the case of the new FO expansions, where the 
approximation is quite good near the spacelike axis and 
gradually deteriorates for points closer to the timelike 



axis. Thus, the new CI expansions considered in the 
present work have a more solid theoretical basis than the 
new FO expansions. 

In Sec. IVIIi we present a determination of as{M^) 
based on the new CI expansions defined in Eqs. ((29)) - 
(1301) . We emphasize that in our analysis we do not rely 
on models and make no assumption about the strength of 
the leading singularities in the Borel plane. The predic- 
tions of the various expansions, reported in (j4ip . exhibit 
a remarkable consistency among each other. 

Our final prediction (j43p , obtained by averaging the in- 
dividual values dH]), is very close to that of the standard 
FO expansion, while the standard CI expansion gives a 
value larger by about 0.022. It is important to under- 
stand the origin of this result. In our opinion, it is re- 
lated to the consistent treatment of the running of the 
coupling and the Adler function coefficients in the stan- 
dard FO expansion and the new CI one. 

As discussed in [2l[, the standard FO expansion is 
suitable for models like the ansatz p3l) - ([55|) . where the 
residues of the dominant renormalons are fixed in a nat- 
ural way from the first coefficients Kn. On the other 
hand, as noticed in Sec. IVIi this expansion is not so effi- 
cient for more artificial models like that presented in Eqs. 
(I57|) - ([M)) . where the strength of the first IR renormalon 
is forced by hand to a lower value. So, the standard FO 
expansion seems more suited than the standard CI one 
for describing functions with a natural pattern of leading 
singularities. Suitable compensations of the two terms 
of the same order in ([3]) playan important role in this 
description. As noticed in [28|, these cancellations are 
destroyed in the standard CI expansion (jlOp , which sums 
the running coupling terms, but drops the Adler function 
coefficients Kn in higher orders. 

It is precisely this deficiency that is corrected by the 
new CI expansions, which sum also the Adler function 
coefficients, by properly implementing the singular be- 
havior near the leading renormalons (with no assump- 
tions about their strength) and expanding in powers of 
a conformal mapping. So, the new CI expansions sum 
both the running coupling terms and the expansion of 
the Adler function, while in the standard FO expansion, 
the fixed-order option is made for both expansions. This 
symmetric treatment explains why their predictions are 
similar, at least for truncation orders = 4 or = 5 
of interest at present. At higher N, while both the stan- 
dard FO and CI expansions start to diverge, the new CI 
expansions show an impressive convergence for all types 
of expanded functions. 

As we already discussed, the largest errors in (|4ip are 
due to the conservative range K5 = 283 ± 283 adopted 
for the 5-loop coefficient K^. Expressed in other words, 
the new expansions appear to exclude both very large 
and very small (or negative) values of K^, which would 
require unusual values of as(M^) to reproduce the input 
(Uni)- For the standard expansions, the error related to 
is much smaller, but it cannot fully account for the 
asymptotic character of the series, which may start to os- 
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cillate at a certain N without approaching the expanded 
function with a sufficient accuracy. 

The value given in (|43| represents our best determina- 
tion, obtained as the average of the determinations (PT|) 
with a very conservative treatment of the uncertainties. 
Our analysis shows that for increasing the precision of 
as(M^) determination with the new expansions a more 
precise knowledge of the 5-loop coeSicient K5 is crucial. 
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Appendix A: Proof of lemmas given in Sec. IIIII 

• Proof of Lemma 1: We remind the reader that 
Lemma 1 was stated and proved in Ref. Q; we recall 
it for completeness here. 

Let us define 



/(Z2)-51(4""(^2)) 



(Al) 



for Z2 G /C2, where z\ ^' is the inverse to Z2, which exists 
since Z2{u) is a conformal mapping. 

The function 7(2:2) is holomorphic on the unit disk /C2 
of the 2:2 -plane and maps this disk into the unit disk /Ci of 
the zi-plane, i.e. \f{z2)\ < 1. Moreover, since zi{Q) = 
and Z2{Q) = by assumption, it follows that /(O) = 0. 

We now apply Schwarz's lemma, which states that if 
a function F{z) is holomorphic on the disk |z| < 1 and 
satisfies the conditions F{0) — and |i^(2)| < 1 for |z| < 
1, then 



\Fiz)\<\z\ 



(A2) 



everywhere in |z| < 1. Besides, if the equality sign occurs 
in (IA2|) at least at one interior point, then it takes place 
everywhere and F{z) has the form F{z) — z exp(za) with 
a real. 

Applying Schwarz's lemma to the function / defined 
in (|A1I) . we have 1/(22)! — kal for Z2 G JC2. Using the 
definition (jAip and the obvious relation 23 ~ u for 

M G 2?2, we obtain 



\Z1{U)\<\Z2{U)1 UGV2. 



(A3) 



Ignoring the mappings that reduce to mere rotations ac- 
cording to Schwarz's lemma, we are left with a sharp 
inequality in (jA3[) . 



|2i(m)| < |22(u)|, ueV2, U^Q, (A4) 
which proves Lemma 1. 



• Proof of Lemma 2: The relations ([18)) imply that the 
coefficients |c„j| can, for large enough n, be represented 
in the form 



— g9j(") 



J 



1,2, 



(A5) 



where gjiji) are real- valued functions, subject to the con- 
ditions liran^oo gj{n) /n = 0, j = 1,2. Then, the ratio 
defined in (fTOl) can be written as 



(A6) 



where 

g{n) = 51(71) - g2{n), 



p(u) = |fi(u)/i2(M)|- (A7) 



Taking the logarithm of (|A6|) . one obtains, for large n, 
the inequality 



InTZniu) 



5(") 



In p{u) 



<0, 



(AS) 



since from (jA7p it follows that lim„_yoo gin-) /^^ = Oj while 
p{u) < 1 for all u G I?2, u ^ Q, according to Lemma 1. 
This implies (jl9p . proving Lemma 2. 
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